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ABSTRACT: The drift and diffusion of a cloud of ions in a fluid are distorted by an inhomogeneous 
electric field. If the electric field carries the center of the distribution in a straight line and the field 
configuration is suitably symmetric, the distortion can be calculated analytically. We examine the 
specific examples of fields with cylindrical and spherical symmetry in detail assuming the ion dis- 
tributions to be of a generally Gaussian form. The effects of differing diffusion coefficients in the 
transverse and longitudinal directions are included. 
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1. Introduction 

We address the time and spatial evolution of the ionization produced along a charged particle's track 
in a fluid (gas or liquid) subjected to a non-uniform, time-independent electric field. The prototypical 
calculation is for a localized ion density produced at some point Ro at time t = 0, described by a delta 
function 5(r — Ro). Of interest is the diffusion in both the direction of the electric field and transverse 
to it, which governs the distributions in arrival time and transverse position at some detecting surface. 
We focus on a volume that is bounded by cylindrical electrodes that produce a radial electric field, but 
begin with more general considerations. Spherical geometry is considered in an appendix. 

The motion of electrons and atomic or molecular ions in electromagnetic fields has a huge ex- 
perimental and theoretical literature. See, for example, the texts and cited literature of Huxley and 
Crompton[[T|] and Mason and McDaniel[^]. An extensive literature about diffusion of charged parti- 
cles in plasmas also exists, but conditions there are far from our concern of relatively few ions moving 
in a neutral medium. 

We make the simplifying assumptions that the ion mobility pL and the longitudinal and transverse 
diffusion constants, Dl and Dj, are independent of electric field strength E or rather of E/N, where 
N is the particle number density of the fluid. The data compilations of transport properties of gaseous 
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ions ^ ||] show that the mobilities of ions in noble gases are field-independent (to ±10 %) typi- 
cally for E/N < 60 — 100 Td.[||] The longitudinal diffusion constants D L for ions in noble gases are 
similarly roughly field-independent typically for E/N < 10 — 20 Td. For this range of E/N, the ratio 
of the diffusion constants is Dl/Dj 1; only at higher values of E/N is Di/Dj greater than unity 
(and both are field-dependent). See, for example, figure 5-2-4 of Ref. [Q]. 

Within our assumptions, we hypothesize a regime in which the microscopic (atomic) processes 
that establish the local drift velocity and diffusion parameters occur on length and time scales small 
compared to the changes in magnitude or direction of the electric field or the average drift position. 
Such a regime permits the use of the simple diffusion equation rather than the full Boltzmann equation. 

For electrons in gases, field-independent mobilities occur typically only for E/N < 10~ 2 Td,with 
various field dependencies at larger £"/A^.[jl]| Diffusion constants for electrons often exhibit rapidly 
varying field-strength dependence because of resonant behaviour in the momentum transfer cross 
sections. Our approximations seem to have limited applicability for electrons. 



2. Preliminaries and initial Ansatz 

Ions in a fluid and subject to an electric field can be described by a charge density Y(r,t) satisfying 
the diffusion equation 

— = ^ + vVii/ = DVV. (2.1) 
dt dt 

We take the velocity v to be given by the mobility /I times the electric field E. We imagine 
that from the prescribed electric field and mobility we can calculate the classical trajectory ro(?) of a 
point charge located at r = Ro at time t = 0. We indicate the velocity along the path by vo(?) and, in 
particular, 

v o (0=iUE(r o (0). (2.2) 

The trajectory is given implicitly as 

r o (0=Ro+ f dt'nE(r {t')). (2.3) 

Jo 

We are seeking a solution proportional to 8(r — Ro) at t = 0. We begin with an Ansatz for y. 

V = [d (4^W2 /2 ex P I" t r ~ *>(0]<[fr + S(t)ij} [r - r (t)]j/(4Dt)} . (2.4) 
The factor [det(/ + g)] 1//2 ensures that y is normalized properly: 

fd 3 ryf=l, (2.5) 
which must be the case if we solve the initial equation, which derives from the continuity equation. In 



the absence of the electric field, the familiar solution is given by Eq.(|2.4[) with gjj = and ro(t) = Ro- 



The basic assumption is that the form is Gaussian with respect to the difference between the point 
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at which the density is evaluated and the location of the classical trajectory at that time, ro(f). The 
form is motivated by the observation that, not only in the absence of the electric field, but also for a 
constant electric field, the solution is given by Eq^^) with gjj = . We calculate the time and spatial 
derivatives of y, writing N(t) for [detfY + g)] 1 / 2 and p, for (r — ro(f)),-: 

4Dt ^- = Pi(S u + gij)pjj-6DY+2p,{d i j + g i j)v jY 

-ptJZpjY+ADt—r, (2-6) 
4DtV k y = -2{8 jk + g jk )pjT, (2-7) 
4D 2 tV k V k y = (8kj + gkj)pj{S H +g k i)pij-2(8 kk + g kk )D\j/ 

= P j (8 j i + 2g jl + g kj g kl )p l ^- -2(3 + gkk)Dr, (2-8) 
4£>f i uE • Vy = -2pLE k {8 kj + g jk )pjY- (2-9) 

We now expand E(r) around ro(t), keeping first order terms, 

v (0 - ME(r) = Ai(E(r o (0) - E(r)) = -|iME(r (0) + • • • • (2-10) 

The omitted terms in the expansion are of order (p /L) n , n>\, times the first order term, where L is the 
length scale over which the electric field changes appreciably. If the spread of the diffusion, 0(p max ), 
is small compared to L, our approximation should be reliable. We combine these into Eq.(2.1) to find 
the equations for gif 

(Sij + gij) , .r^lnN dgij 

Pi^-^PJ - eD + ADt—y-p^pj 

= pj (S jI + 2 gj , + g kjgkl ) pi _ 2{3 + gkk)D 

+2np i d i E k (S k j + gj k )p j . (2.11) 

The terms with and without quadratic dependence on p must cancel separately. In fact, we find that 
the equations that derive from terms with quadratic dependence ensure that the terms without that 
dependence cancel, as well. The terms without the quadratic dependence give 

It d\nN 1 
l- T -r— = 1 + ^ (2-12) 
3 dt 3 



The terms with quadratic dependence on p give 

(Sjj + gjj) dgjj _ (5j i + 2gj i + g kj g ki ) 
t dt t 

or 



+ 2nd i E k (8 kj + g jk ) (2.13) 



dy_ + gij + 8*gkj +2 d . m + ) = (2 . 14) 
dt t 
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3. Restriction to rectilinear motion and symmetric geometry 



Only the symmetric part of gjj enters \jf. Consistency requires that the antisymmetric part of Eq.( 2.14 ) 
vanish: 

(Sij + gij)d k Ej = (8 kj + gkj )diEj (3. 1) 

The vanishing of the curl of E is enough to take care of the 8{j portion but we are left with the 
requirement 

gijdkEj = gkjdiEj (3.2) 

In general, the product of two symmetric matrices is not symmetric. We need to impose additional 
constraints. In particular, we suppose that the motion is rectilinear along the x direction, i.e. r(t) = 
r(t)x, E(?) = E x (t)x. Write the electrostatic potential in the neighborhood of the line of motion as an 
expansion in y and z: 



®(x,y,z) = <l>(x,0,0)+y&y(x)+z& z (x) 
+\y 2 <$>yy{x)+ l -z 2 <$> zz {x) 
+xy®xy(x) +xz&xz(x) +yz%z{x) + ■ ■ 



(3.3) 

By hypothesis the electric field along y = 0, z = is entirely in the x direction. Then ^(jc) = 



& z (x) = &xy{x) = *&xz(x) = an d Eq.(3_3) becomes, to second order inclusive, 



1 7 1 7 

&(x,y,z) = 4»(x,0,0) + -y 2 ®yy {x) + -z 2 <t> zz {x) +yz& yt {x) 



(3.4) 



If the surface charge distributions determining the electric field are symmetric under y — ► — y and 
z — > — z, the term yz& yz vanishes along with <t> ) ,<t> z ,<t> X} , and <I> VZ , and the only non-zero values of 
d,Ej along the line (x, 0, 0) occur for i = j. Two clear examples are cylindrical symmetry with the z, 
direction taken as the axis of the cylinder and spherical symmetry. In the former case all derivatives 
with respect to z vanish so d z E z = and in the latter, the derivatives with respect to y and z are equal 
so d z E z = d y Ey. See figure [l| 

Under any of these circumstances we can write 

d i E j = C i (x)8 ij (3.5) 



and we can satisfy Eq.(|3.2|) with the ansatz 

gij = gi$ij (3-6) 



i.e. both gjj and djEj are diagonal along the linear trajectory path of Eq.(2_3). Now Eq.( 2.14 ) simplifies 
to 

djjt + Si^+Si) + 2iidiEi{\ +gi) = (3.7) 
where there is no summation on i. 
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Figure 1. The Ansatz for Eq.(2.4), is applicable for rectilinear trajectories if the matrix djEj is diagonal, 
for a suitable choice of coordinates. Taking x along the direction of motion, a charge distribution (or potential) 
symmetric under y ~ * —y and z — » — z will ensure this. The figure shows such a trajectory and charge distribu- 
tion. Here only the trajectory shown has the required symmetry. For the cylindrical and spherical geometries 
all radial trajectories have the symmetry required by our Ansatz. 



4. Example of cylindrical geometry 

For a cylindrical geometry, uniform in the z direction and having azimuthal symmetry, with a radial 
electric field, we have 

E = -f, (4.1) 
r 

where r = xx + yy is now two-dimensional. We assume here that the signs of the charge and field are 
such that the force is radially outward. 

We can anticipate the effect of the inhomogeneous electric field, Eq.(fi~l|). The portion of the 
cloud that is further out than the center of the cloud experiences a weaker field and thus moves more 
slowly than the outwardly moving center, while the portion of the cloud closer in to the central axis 
experiences a greater field and moves more rapidly than moving center. Thus the spread of the cloud 
in the radial direction should be smaller than in the absence of the electric field. On the other hand, 
the diverging electric field will tend to spread the cloud transversely, leading to a larger dispersion 
than in the absence of the electric field. 

Take the initial position Ro = (#0,0,0). The motion is along the x axis where 

A A 

d x E x = — d y Ey = ; d x E y = = d y E x , (4.2) 
x L yr 

satisfying the symmetry requirements for our Ansatz. Thus we have the diagonal form djEj = Cj(x)8jj 
where 

Ci(jc) = ~ (4.3) 
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C 2 (x) = A (4.4) 

C 3 (jc) = 0. (4.5) 

The trajectory is given by ro(t) = ro(t)x, with 

r (t) 2 =xl + 2\xAt 1 (4.6) 

so that 



2jUd A £ v = - ? = -2]xd y E y . (4.7) 



2jUA 
^ + 2juAf 

Introduce 

T = A (4.8) 

so that 



ro(t)=xoJl + ^. (4.9) 

In the characteristic time T the diffusion center would go a distance jco/2 if the field had the constant 
value E =A/xq. When t = 3T, the classical particle is twice as far from the axis of the cylinder in the 
l/r field as it was at t = 0. From Eq.(3.7): 



Finally, writing t = Ts 



dgi gi(l+gi) , 1 ( . , v 
at t T + t 

dgi g2(l+g2) 1 



dt t T + t 

dg3 g3(l+g3) 
dt t 



The equation for #3 has the solution 



;i+g2); 



dgi gi(l + gi) , 1 \ . 

dgi g2(l+g2) 1 



ds s \+s 

dg3 _ g3(l+g3) 
ds s 



(i+g 2 ); 



(4.10) 



(4.11) 



g3 = -^-r (4.12) 
Bs — 1 

where B is a constant. When s = 0, all the gt are zero. It follows that B = 00 and #3 remains zero: 
Diffusion in the z direction is unchanged by the presence of the radial electric field. 
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Figure 2. The functions appearing in the Gaussian description of the charge density y/ for outward drift, as a 
function of the dimensionless time variable s = t/T for cylindrical geometry. 



Eqs.( ft.ll| ) for g\ and g2 are Ricatti equations and can be solved analytically. The details of the 
solutions are given in Appendix A. The results are 



8\ 



2 + s 
82 = -1 + 



(l+j)ln(l+j) 



(4.13) 
(4.14) 



The behavior of the functions gi and g2 is shown in figure |2[ 

The effective parallel and transverse spreads relative to On = V^Dt are 



1 Oj_ 



1 



shown in figure ||. In the z-direction the spread is 03 = Co- 



(4.15) 
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Figure 3. The effective dimensionless widths Eq.(4.15) for outward drift, in the directions parallel and perpen- 
dicular to the radial direction as functions of the dimensionless time variable s — tjT for cylindrical geometry. 
If the diffusion constants differ in the longitudinal (parallel to the field) and transverse directions, multiply the 
ratio o]|/ob by \/Dl/Dt and take Ob = yJIDjt. (see Sec. 6). 
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Figure 4. The functions appearing in the Gaussian description of \jf for inward drift in the cylinder. Negative 
values of s correspond to positive values of t according to t = Ts. For s = — 1, ro(t) = 0. 

5. Inward drift 

We have imagined that the radial drift is outward, towards larger r. If the sign of the electric field is 
reversed, the derivation goes through much the same way except that the sign of T is negative. Since 
t = Ts, positive times then correspond to negative values of s. The center point of the diffusion cloud 
arrives at the central axis at s = —1. Because the field increases in the direction of motion, we expect 
the effect of the field to be opposite to what we found above: the dispersion should be greater in the 
radial direction and smaller in the transverse direction. 

We show the results obtained for inward drift in Figs. |] and H. 

6. Different coefficients of diffusion for longitudinal and transverse directions 

It is straightforward to include the effect of having different diffusion coefficients Dl and Dj for the 
longitudinal and transverse directions. For the Ansatz we take 




(6.1) 
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Figure 5. The effective dimensionless widths Eq.( 4.15 ) in the directions parallel and perpendicular to the radial 
direction for inward drift in the cylinder. Negative values of s correspond to positive values of t according 
to t = Ts. If the diffusion constants differ in the longitudinal (parallel to the field) and transverse directions, 
multiply the ratio On / Ob by \[~DiJDt and take Ob = (see Sec. 6). 



with Di = Dl,D2 = Dt, = Dj. Exactly the same differential equations emerge, but setting ah = \J2Djt 
we have 




1 



31 l°L i 
On 

Oj_ 

Co y/l +g2 

2 = 1 



(6.2) 
(6.3) 
(6.4) 



7. Diffusion of a track segment 

To illustrate the possible use of our results, instead of a single point source arising at r = Ro at t = 0, 
we consider the ionization along the path of a relativistic charged particle through a cylindrical drift 
chamber. The ionization trail is assumed linear (in both senses) and is, on the time scale relevant 
for the subsequent diffusion, deposited instantaneously. The resulting charge distribution is thus a 
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Figure 6. Left (a): Track segment (Ri,R 2 ) with typical starting point R(4). Right (b): Details of a typical 
starting point showing its decomposition R(§ ) = Rx (%)+R\\z into a piece along the drift direction and a piece 
parallel to the axis of the cylinder. 



superposition of distributions of the form Eq.(JO|) for contributions along the particle's track. For 
definiteness we consider a track segment with end points Ri, R2, as shown in figure (|6]a), with a 
typical point on the track labeled by R(£). For convenience, define the distance between the end 
points as L = |R 2 — Ri|. Then we have R(£) = Ri + (E, /L)(R 2 - Ri). In figure (§b), the typical 
starting point R(£) is singled out. We define R\\{E,) = z-R(£) and R±(^) = R(£) -zR\\{£,). Then 
the drift center's coordinate r (f) can be written according to Eq.(^) as 

r (f) = R^l+t/T^)] 1 / 2 + zR { ^) (7.1) 

where from (30), since xq(^) = |Rj_(£)|, we have 

r«) = « (7-2) 

In the exponent of \j/ the Cartesian components of the coordinate difference, (r — ro(t)), must be taken 
with respect to Rj_ ) as the x-direction: 

Pl = (r-r o (0)i = (r-ro(O)-ftx(O (7.3) 
p 2 = (r-r (0) 2 = (r-r (())-(zxRi(O) (7.4) 
P3 = (r-r (0)3 =z-R\\{£) (7-5) 



The normalized diffused charge density in space and time for the track segment can thus be 



written as a numerical integral over y given by Eq.(|6JJ), 



¥(r,f|Ri,R 2 ) 



(7.6) 
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The functions gi(s) and g 2 (s) are to be evaluated at s = t/T(£); g3<V) = 0. 
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A. Solution of Ricatti equations for cylindrical geometry 

The procedure for solution of the Riccati equations can be found in Refs.[||, [J. The general form is 

y = qo(x)+qi(x)y + q 2 (x)y 2 . (A.l) 
To find the solution we consider the associated second order homogeneous linear differential equation 

u" -P(x)u +Q(x)u = (A.2) 

with 

Q = q 2 q Q ; P = q\ + 

<?2 



Q = q 2 qo\ P = qi + ^. (A.3) 



Then the solution to the original equation is 



u' 

y = ■ (A.4) 

q 2 u 

(a) For g\ we have 

1 1 1 

qi = ; qi = — 7z r ; qo = "] ; > ( A - 5 ) 

S S(l+ S) 1+5 

1 2 + s 

^(1+5) 2(1+5) 

and 

„ 2+s , 1 
u + — -u tU = 0, 

5(1+5) 5(1+5) 

s{\+s)u" + {2 + s)u -u = 0. (A.7) 
We recognize this as the hypergeometric equation for a function w(z)[Q], 



d 2 w , , . ,dw 

— + [c-(a + b + m- 

whose solution is 



z(l-z)— + [c-{a + b+\)z}— -abw = 0, (A.8) 



, , ab a(a + l)b(b + 1) z 2 

c c(c+l) 2! 

In particular, we have 

Mi =F(1,-1;2;-j) = 1 + -5 (A.10) 
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and finally, from Eq.( |A.4D , we have 
(b) For g2 we have 



12 



so 



and 



1 



Q 



2 + s 
2s+l 



1 



s(i+ s y i+.v" 

1 „ 2 + 3s 



s(l+s)' 2(1+*)' 
„ 2 + 3s , 1 



(A. 11) 

(A. 12) 
(A.13) 



a + 



-u + ■ 



-u = 0, 



3(1+3) s(l+s) 
s(l+s)u" + (2 + 3s)u! + u = 0. 



with the solution 



u 2 = F(l,l;2;-s) = --ln(l+s), 



from which we find, from Eq.( |A.4| ) 



g2 = -l + - 



(A. 14) 
(A. 15) 

(A. 16) 



(l+j)ln(l+j) 

This can be confirmed by direct substitution in the original differential equation. 
B. Example of spherical geometry 

The volume consists of the region between two concentric conducting spherical shells of radii a, b (a< 
b). Voltages are applied to the shells such that a radial electric field E = Ar/r 3 is produced. It is as- 
sumed that the force on the ions is outward, unless stated to the contrary. 

Consider the starting point to be at R = (R, 0, 0) at t = 0, with a < R < b. The drifting center moves 
outward along the x =x\ axis. The classical equation of motion, v = pL E, reduces to dx/dt = [lA/x 2 , 
with solution, 

r (0 = {R 3 + 3nAt) 1/3 



We define 



Then we have 



*(1 + f) 1/3 



(B.l) 
(B.2) 
(B.3) 



The needed partial derivatives of 2/1 E near the x -axis are 



2\x d\E\ = - — 
2\x d 2 E 2 = + 



1 



3T (1 + t/T) 
1 1 



3T (1 + t/T) 
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with the derivative in the #3 direction the same as for x%. The equations for gj, j = 1, 2 take the form 

d_sj + Sj (i+ g j) +a (A±iA = (B . 4) 

where, as before, s = t/T. The parameter a is a = —2/3, +1/3 for the longitudinal (7 = 1) and 
transverse (j = 2) directions, respectively. 

With the introduction of gj = s- w(s), the equation for gj becomes a Ricatti equation for w(s): 

dw 7 .2 a . a „ „ _ x 

^ + W + w( 7 + T+^ ) + ^y = a (B - 5) 



With z. = — s and w = this becomes the hypergeometric equation, Eq^AJS), for F(a,b;c;z) with 



c = 2, a + Z? = 1 + a, ab = a. Thus for j = 1, a = 1, b = —2/3, c = 2 and for j = 2, a = 1, ft = 
1/3, c = 2. With some algebra and use of relations of the hypergeometric functions, we find the 
solutions for the longitudinal and transverse functions, g\ and g2 : 

_ 1 ; F(l, 1/3; 3; */(!+*)) 
gl 3 (l+s) F(l; -2/3; 2; */(! + ■*)) 1 ' 

Zl g F(l,4/3; 3; s/(l+s)) 
82 6 (l+s) F(l; 1/3; 2; */(! + *)) ' 
The classical center is given by ro(t) = R(l + s) 1 / 3 . For example, for ro(t) /R = 2, we have s = 7. 

Figures ^ and [8] show the behavior of gi(s) and — g2(s) versus s and the corresponding widths, 
Oli/oo and <7j_/ob> qualitatively similar to the cylindrical situation, but with differences stemming 
from the more rapid fall-off of the electric field with distance. 
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Figure 7. The functions appearing in the Gaussian description of \jf for spherical geometry and outward drift 
as functions of the dimensionless time variable s = t/T, where T is now given by Eq.(B.2). 
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Figure 8. The effective dimensionless widths in the directions parallel and perpendicular to the radial direction 
for spherical geometry and outward drift as functions of the dimensionless time variable s = t/T, where T is 
now given by Eq.(B.2). 
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